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INRIA49.FRD 

nispositif de simulation de deformations — de — mater iaux, 
notamment de tissus co rporels mous 

L* invention concerne la simulation de deformations de 
materiaux, notamment de tissus corporels mous. 

II est prevu habituellement une interface de visualisation 
pour representer au moins un objet suivant une forme predeter- 
mine, maillee avec un motif de maille choisi. Une interface 
utilisateur munie d'un organe mecanique de saisie permet de 
simuler une ou plusieurs forces exercees globalement sur 
1' objet. 

Une zone-memoire est prevue aussi pour stocker : 

- des parametres mecaniques du materiau dans lequel 1' objet 
consider est prevu, 

- des donnees de position de 1' objet, prises en des sommets 
d'un maillage choisi, et 

- des donnees de force, representant , en intensity et 
position, une sollicitation a exercer sur 1' objet. 

Un calculateur capable de cooperer avec la zone-memoire eyalue 
de nouvelles positions des sommets, en fonction de la 
sollicitation exercee et des parametres mecaniques du 
materiau. Pour visualiser la deformation et/ou le retour a une 
position de repos de 1' objet en mouvement, le calculateur 
determine une evolution temporelle des positions des sommets 
en fonction de cette sollicitation. Dans une application a la 
simulation de la deformation (ou de la decoupe) de tissus 
corporels, un personnel qualifie (tel qu'un chirurgien) peut 
suivre la simulation d'une deformation et/ou un retour a une 
position de repos . 

Cependant, les deformations simulees avec de grands deplace- 
ments ne sont pas fiables habituellement. 

La presente invention vient am61iorer la situation. 
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Selon une premiere caracteristique important de 1' invention, 
le calculateur comprend: 

- un module apte a determiner iterativement , pour chaque 
maille, l'ecart entre la longueur courante d'une arete de la 
maille et sa longueur au repos, ainsi qu'a stocker des donnees 
de force respectives pour chaque sommet de la maille, et 

- un module apte a determiner iterativement, pour chaque 
sommet, de nouvelles donnees de position de ce sommet en 
fonction de la composition des forces qui s'y exerce, de 
donnees portant sur au moins une position anterieure du sommet 
et des parametres mecaniques du mater iau. 

Dans une forme de realisation preferee, le motif de maillage 
choisi comprend au moins une base triangulaire. 

Ainsi, si l'objet presente la forme d'une enveloppe tridimen- 
sionnelle creuse, le maillage choisi est pref erentiellement 
triangulaire, de sorte que le calculateur est apte a determi- 
ner une composition des forces en chaque sommet de triangle, 
en fonction de l'ecart entre la longueur courante de chaque 
arete du triangle et la longueur de cette arete au repos. 



Si l'objet presente une forme tridimensionnelle pleine, le 
maillage choisi est pref erentiellement tetraedrique, de sorte 
25 que le calculateur determine une composition des forces en 
chaque sommet de tetraedre, en fonction de l'ecart entre la 
longueur courante de" chaque arete du tetraedre et la longueur 
de cette arete au repos. 
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Avantageusement, la zone-memoire est agencee pour stocker, en 
association avec chaque maille, des parametres mecaniques du 
materiau de la maille, au moins partiellement definis 
localement, en particulier au niveau de la maille ou d' ele- 
ments de celle-ci. 

Dans la demande publiee FR-9714506 du Demandeur, la simulation 
des deformations de l'objet etait, dans la realisation 
decrite, basee sur un modele lineaire. Bien que 1' utilisation 
du dispositif decrit dans FR-9714506 s'est averee tres 
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prometteuse dans le domaine de 1' assistance chirurgicale pour 
la simulation de la decoupe de tissus corporels mous (foie, 
peau, etc), la simulation de grands deplacements , le cas 
echeant non lineaires, n'etait pas completement naturelle. De 
plus, ce dispositif ne permettait pas de prevoir la reaction 
a une force exercee sur un type quelconque de materiau 
deformable, en particulier pour de grands deplacements. 

Dans une realisation preferee de la presente invention, le 
calculateur est agence pour determiner ledit ecart de 
longueurs courante et au repos, en vue d'estimer une derivee 
de l'energie potentielle de deformation de chaque maille, 
exprimee en fonction d'un tenseur de Green-St Venant et de 
coefficients mecaniques propres au materiau dans cette maille. 

La derivee de l'energie potentielle, ainsi exprimee, par 
rapport a la position du sommet, permet d'obtenir la force (ou 
la composition des forces) qui s'y exerce. 



20 



Dans 1' estimation des derivees de l'energie potentielle aux 
sommets, exprimee ainsi en fonction d'un tenseur -de Green- 
St Venant, le calculateur est agence pour determiner, plus 
precisement, des differences entre les carres des longueurs 
courante et au repos de chaque arete, en vue de determiner la 
25 composition des forces precit6es. 

Les parametres mecaniques precites comprennent avantageusement 
les coefficients de Lame du materiau dans chaque maille 
consid6r£e . 
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Le dispositif selon 1' invention permet alors de simuler des 
grands deplacements et/ou d'^lasticite non lineaire, le cas 
echeant d'objets heterogenes, realises dans des materiaux 
dif f erents . 

II est possible ainsi de simuler la deformation d'un objet 
tridimensionnel plein, le cas echeant enveloppe dans une peau 
creuse. II conviendra en l'espece d'appliquer un maillage 
tetraSdrique a 1' objet et un maillage triangulaire a la peau 
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et d'affecter a chaque type de maille des parametres mecani- 
ques adequats. H peut etre prevu aussi de simuler les 
deplacements d'objets mous (tels que des muscles), en presence 
d'objets rigides et solidaires (tels que des os), comme on le 
5 verra plus loin. 

La presente invention peut se manifester aussi sous la forme 
d'un procede de simulation de deformations d'objets, notamment 
de tissus corporels mous, comprenant tout ou partie des etapes 

10 suivantes : 

a) prevoir une forme predetermine de l'objet au repos, 

b) mailler ladite forme, sur la base d'un motif choisx, 
preferentiellement de base triangulaire (triangle ou te- 
traedre ) , 

15 c) memoriser les donnees de position des sommets des maxlles 
dans l'objet au repos, 

d) simuler une ou plusieurs forces exercees globalement sur 
l'objet (ou des deplacements de l'objet), 

e) estimer, en fonction des proprietes mecaniques du mater iau 
20 et des positions des sommets precedentes, une force resultante 

sur chaque sommet de chaque maille, et 

f) en fonction de la force exercee sur chaque sommet, 
determiner le deplacement et, de la, la deformation et/ou le 
retour a une position d'equilibre des sommets des mailles. 

25 

L'etape c) passe en outre par une attribution d' indices des 
elements des mailles, dans le cadre de 1' elaboration d'une 
structure de donnees, et, le cas echeant, par une simplifica- 
tion de la structure de l'objet, comme on le verra plus loin. 
30 Les etapes e) et f ) s'appuient sur un modele developpe par le 
Demandeur et que l'on decrira en detail ci-apres. 

D'autres caracteristiques et avantages de 1' invention 
apparaltront a l'examen de la description detaillee ci-apres, 
35 et des dessins annexes sur lesquels : 



- la figure 1 represente schematiquement un dispositif 
une forme de realisation de 1' invention, 
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- les figures 2A et 2B representent schematiquement des 
mailles triangulaires, respectivement initiale et deformee, 
avec les notations des points du triangle utilisees par la 
suite, dans 1 ' estimation de la composition des forces en 
chaque sommet du triangle, 

- la figure 3A est un schema de tetraedre aidant a comprendre 
les notations utilises pour un point Pj du tetraedre et le 
vecteur normal nj de la face oppos^e, 
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- les figures 3B a 3D sont des schemas de tetraedres aidant 
a comprendre, avec la figure 3A precitee, comment se calculent 
des sommes sur des aretes et des couples d' aretes d'un 
tetraedre, pour ^valuer une force que subit un sommet Pp du 

15 tetraedre, 

- les figures 4A et 4B sont des schemas d'aide a la comprehen- 
sion de 1' attribution d' indices dans differentes mailles 
triangulaires , 

20 

- les figures 5A et 5B sont respectivement des schemas de 
motifs triangulaires comprenant un sommet courant (repr6sente 
par un carre) et un sommet virtuel ( represents par un disque), 
avant simplification, et de la meme forme maillee apres 

25 simplification, 

- la figure 6 represente schematiquement un exemple d'un 
maillage triangulaire non conforme, et 

30 - la figure 7 represente un organigramme sch^matique du 
deroulement des calculs de forces et de positions des sommets 
represent ant l'objet. 

La presente description fera mention de formules et d'equa- 
35 tions qui sont retranscrites en annexe. . 

Les dessins, ainsi que 1' annexe, contiennent, pour l'essen- 
tiel, des elements de caractere certain. lis pourront done non 
seulement servir a mieux faire comprendre la description, mais 
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aussi contribuer a la definition de 1- invention, le cas 
ech6ant. 

On se refere tout d'abord a la figure 1, sur laquelle le 
5 dispositif se presente, dans la realisation decrite, sous la 
forme d'un ordinateur PC. L'ordinateur comprend une memoire 
permanente MEM (par exemple de type ROM) dans laquelle sont 
enregistres des programmes de calcul des positions des sommets 
de maille et les forces resultantes en ces sommets, ainsi 
10 qu'un processeur uP capable de cooperer avec une memoire de 
travail MT (par exemple de type RAM) pour effectuer ces 
calculs* 

Les modules precites du calculateur decrit plus haut peuvent, 
15 le cas echeant, se presenter sous la forme de programmes 
stocked dans la memoire MEM, ou encore etre charges dans la 
memoire de travail MT, a partir de la lecture d'un support 
dans un peripherique externe (non represents ), tel qu'un 
lecteur de disquettes, de CD-ROM, ou encore recus a partir 
20 d'une liaison de communication, par exemple par MODEM. Bien 
entendu, de tels programmes apparaissent comme des moyens 
particulierement importants pour mettre en oeuvre 1' invention. 
A ce titre, la presente invention vise aussi de tels program- 
mes, le cas echeant sur un support du type precite (disquette, 
25 CD-ROM, etc), dont des exemples de realisation detailles 
seront donnas plus loin. 

L'ordinateur comprend en outre un moniteur VISU relie au 
processeur uP par 1 ' intermediate d'une interface graphique 

30 IG, 1' ensemble VISU et IG (appele "interface de visualisation" 
ci-avant) permettant de visualiser la forme d'un objet dont 
la deformation est simulee. L'ordinateur comprend encore une 
interface utilisateur IU relie a un organe mecanique de saisie 
tel qu'un clavier CLA ou une sour is MO, permettant a un 

35 utilisateur de simuler une force globalement exercee sur 
1' ob jet, ou encore un deplacement, avec deformation initiale, 
de cet objet. 
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aussi contribuer 2i 1 ' interpretation de l f invention f le cas 
echeant . 



On se ref^re tout d'abord 2t la figure l r sur laquelle le 
5 dispositif se prSsente, dans la realisation d6crite, sous la 
forme d'un ordinateur PC* L ' ordinateur comprend une memoire 
permanente HEM (par exemple de type ROM) dans laquelle sont 
enregistres des programmes de calcul des positions des sommets 
de maille et les forces resultantes en ces sommets, ainsi 
10 qu'un processeur uP capable de coop^rer avec une m6moire de 
travail MT (par exemple de type RAM) pour effectuer ces 
calculs. 

Les modules pr6cit£s du calculateur deer it plus haut peuvent, 
15 le cas echeant, se presenter sous la forme de programmes 
stock6s dans la memoire MEM, ou encore etre charges dans la 
memoire de travail MT, ci partir de la lecture d'un support 
dans un p6riph6rique externe (non represents ), tel qu'un 
lecteur de disquettes, de CD-ROM, ou encore regus h partir 
20 d'une liaison de communication, par exemple par MODEM* Bien 
entendu, de tels programmes apparaissent comme des moyens 
particuli£rement importants pour mettre en oeuvre 1' invention. 
A ce titre, la presente invention vise aussi de tels program- 
mes, le cas echeant sur un support du type precite (disquette, 
25 CD-ROM, etc), dont des exemples de realisation detailies 
seront donnas plus loin. 

L ' ordinateur comprend en outre un moniteur VISU relie au 
processeur uP par 1 ' intermediaire d'une interface graphique 

30 IG, 1' ensemble VISU et IG (appeie "interface de visualisation" 
ci-avant) permettant de visualiser la forme d'un ob jet dont 
la deformation est simuiee. L' ordinateur comprend encore une 
interface utilisateur IU relie h un organe mecanique de saisie 
tel qu'un clavier CLA ou une souris MO, permettant k un 

35 utilisateur de simuler une force globalement exercee sur 
l'objet, ou encore un deplacement, avec deformation initiale, 
de cet objet. 
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La forme de l'objet a mailler peut etre une surface repre- 
sentee en trois dimensions. Dans ce cas, le maxllage est 
triangulaire (figure 2A) . La forme de l'objet a mailler peut 
etre une tridimensionnelle pleine et dans ce cas, le maillage 
est tetraedrique (figure 2B) . Le cas echeant, un module de 
maillage est memorise dans la memoire MEM et peut etre mxs en 
oeuvre par le processeur uP, en cooperant avec 1' interface 
graphique IG, pour visualiser une forme de l'objet maillee sur 
le moniteur VI SU. 



Ce qui suit decrit un modele developpe par le Demandeur et 
applique a un dispositif de simulation de deformation. Dans 
ce modele, on tient compte d'une forme de repos, sans 
raisonner par rapport a une position de repos (raisonnement 
5 en deplacements qui est habituellement la base des modeles 
classiques) . 

On considere d'abord une surface representee en trois 
dimension, maillee selon un motif triangulaire. Un triangle 
20 TO = (P0,P1,P2) de la forme au repos, d'aretes (l 0 ,li,l 2 )' Ce 
triangle (de gauche sur la figure 2A) subit une transformation 
O pour donner le triangle T = (Q0,Q1,Q2) (a droite de la 
figure 2B). 

25 L'energie potentielle elastique de cette deformation est 
exprimee ici a partir du tenseur de deformation de Cauchy- 
Green C (formule Fl de 1' annexe). En considerant 1' element 
fini obtenu en associant a ce triangle un ensemble de 
fonctions d' interpolation lineaires, le gradient de la 
30 deformation s'exprime en tout point X du triangle par la 
formule F2, ou A est l'aire du triangle TO, nj est le vecteur 
unitaire perpendiculaire au segment (Pj+l,Pj+2), et les aj 
sont les vecteur s de forme du triangle. 

35 Le tenseur de deformation C devient alors donne par 1' ex- 
pression de la formule F3. Avantageusement , ce tenseur de 
deformation est invariant par translation et ne depend done 
pas du choix de l'origine du repere. Preferentiellement , 
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l'origine du repfere est placee au centre du cercle circonscrit 
au triangle deforme T, de rayon r (figure 2A) f ce qui permet 
d'exprimer le tenseur C sous la forme donnee §t la formule F4 f 
ou Lp est la longueur de 1/ arete defonruSe joignant les sommets 
> Qj et Qk. 

L f expression de l*energie potentielle de deformation (donnee 
par la formule F5 de 1' annexe) utilise celles de la trace du 
tenseur C (donnee par la formule F6) et de la trace du tenseur 
0 C 2 (que l'on decrira plus loin). 

k partir des relations classiques dans le triangle TO (loi des 
cosinus), qui s'expriment par la formule F7, on obtient les 
produits scalaires entre les vecteurs directeurs, donnes par 
5 la relation Rl de 1' annexe. Les traces respectives des 
tenseurs C et E s'6crivent alors selon les formules 8 et 9 
donn#es en annexe. Ce terme ne depend que des variations de 
longueur des ar§tes du triangle et de sa forme au repos 
(representee par les grandeurs Ki). 

De l' expression de tr C 2 (formule F10), on deduit trE 2 et, 
finalement, l'^nergie potentielle eiastique de deformation du 
triangle est donnee par la relation R2 de 1' annexe. 

25 La force eiastique appliquee par le triangle sur un de ses 
sommets s' obtient en derivant l'energie eiastique par rapport 
a la position de ce sommet. Dans l'exemple represente sur la 
figure 2A, le sommet PO a pour position deformee QO. En 
derivant la relation R2 par rapport h. QO , on commence 

30 pref6rentiellement par deriver les termes L 2 p pour obtenir 
1' expression de la force suivant la formule Fll f qui est 
1» equation de la force eiastique exercee par le triangle sur 
chacun de ses sommets. Cette force s'exprime bien uniquement 
en fonction des variations de longueur des aretes du triangle 

35 (L^ - li 2 ), des paramfetres mecaniques du materiau (X et yi) et 
de la forme du triangle au repos (donnee par les termes Ki). 



1er depot 

it 

9 

Avantageusement , cette force contient implicitement une 
propriete d'invariance par rotation puisque tous ses termes 
restent nuls tant que les longueurs des ar§tes du triangle ne 
changent pas, c'est-Si-dire tant que le triangle ne se deforme 
5 pas- Ainsi, les calculs de deplacements peuvent se simplifier 
du fait qu'il n'est plus necessaire de compenser les transfor- 
mations rigides, 1'invariance par transformation rigide etant 
intrinsfeque k la formulation du modfele deformable. 

10 En etendant ce qui precede h une forme tridimensionnelle 
pleine (cas volumique ) , maillee avec des tetraedres, la 
formulation de l'#lasticite selon Saint Venant-Kirchhof f 
permet d'exprimer le gradient de la deformee (formule F12), 
ainsi que le tenseur de deformation de Cauchy-Green (formule 

15 F13), ou Aj est l'aire de la face oppos6e au sommet P j , nj est 
la normale a cette face et V est le volume du tetrafedre 
(figure 3A) . 

En plagant le centre du repere au centre de la sphere 
20 circonscrite au tetrafedre T, on peut exprimer le produit 
scalaire des vecteurs position de deux sommets du tetra^dre 
en fonction de la longueur de l'ar§te qui les joint et du 
rayon r de la sph&re, ce qui permet d*ecrire le tenseur de 
deformation selon la formule F14. 

25 

L'energie elastique de deformation s'exprime en fonction du 
tenseur de deformation de Green-Saint Venant (reference E), 
mais, dans une maille 3D (formule F15). L'expression de la 
trace du tenseur C s'exprime encore avec les produits 
30 scalaires (aj.ak), en utilisant la loi des cosinus pour les 
angles diddres et en fonction des longueurs des arStes du te- 
trafedre (formules 16 et 17). Aussi, la trace de C 2 s'ecrit 
selon la formule F18. Finalement, le Demandeur montre que 
I'energie elastique de deformation s'exprime par la formule 
35 F19 de 1' annexe (compietee par la formule F20). 

L 1 expression de la force elastique exercee par le tetraedre 
sur un de ses sommets s'obtient en • derivant I'energie 
elastique par rapport Si la position de ce sommet. Une 
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formulation facilement utilisable s'exprime par la formule F21 
de 1' annexe et s'interprete a l'aide des figures 3A a 3D. 

Dans la formule F21 : 

- dans la deuxieme ligne, k et m sont les indices complemen- 
taires de j et. p, 

- dans la derniere ligne, m est l'indice complementaire de 
j, k et p, et : K jtanp = K jp K^h- K jm K kp - K jk K,^ 

La somme sur les aretes (j,p) signifie qu'il faut parcourir 
les trois aretes adjacentes au sommet Pp (sur lequel s'exerce 
la force) tandis que la somme sur les aretes (j ,k) (indices 
differents de p) fait intervenir les trois aretes situees sur 
la face opposee au sommet Pp. 



En se referant aux figures 3A a 3D : 

- la somme sur les aretes (j,p) fait intervenir trois termes 
(figure 3A) , 

- la somriie sur les couples d' aretes ( ( j ,P) r (k»P) ) fait 
!0 intervenir trois termes (figure 3B), 

- la somme sur les couples d' aretes ( ( j ,p) , ( j »k) ) fait 
intervenir six termes (figure 3C) et 

- la somme sur les couples d' aretes ( ( j ,k) , (m,p;) ) fait 
intervenir trois termes (figure 3D). 

25 

II apparait encore une invariance par rotation des expressions 
ci-avant, pour une forme tridimensionnelle pleine, ce qui 
permet avantageusement d'eviter 1' apparition de distorsions 
irirealistes lorsqu'une partie de l'objet subit une rotation 

30 trop importante. Cette invariance s'etend aussi a la transla- 
tion. Ainsi, l'objet peut subir un mouvement global que lconque 
de translation et de rotation tout en conservant une energie 
et des contraintes internes nulles. Le Demandeur a verifie 
que, dans ce cas, la composante rigide de la deformation est 

35 nulle. 



De plus, la discretisation de l'objet en utilisant des 
elements finis revient a approximer 1* energie elastique sur 
un ensemble de sous-domaines definis par des formes g6om6tri- 
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ques simples* Le module ci-avant propose ainsi, sur chaque 
sous-domaine, une formulation de l'energie qui ne depend que 
du changement de forme de ce domaine. 

5 Avantageusement, la visualisation d'une deformation simuiee 
en considerant une forme de repos initiale (plutdt que des 
positions de repos initiales) parait plus naturelle que celles 
simuiees avec les modules connus. De plus, la consideration 
d'une forme de repos initiale evite les problemes de precision 

.0 qui apparaissent lorsque le modele s'eioigne notablement de 
sa position de repos* 

Ce qui suit decrit des preiiminaires ^ une mise en oeuvre des 
modules pr^cites du calculateur pour la determination des 
15 forces en chaque sommet du maillage choisi. 

II convient tout d'abord de definir une structure de donnees. 
Le Demandeur s'est d'abord confronte §l la difficult^ d'attri- 
buer des indices adequats aux elements de chaque maille- 

20 

Pour des motifs triangulaires , la structure de donnee 
correspondant h une triangulation consiste en une liste de 
triangle , une liste d' aretes et une liste de sommets : 

- pour chaque sommet, on determine une arete adjacente, 

25 - pour chaque arete, on determine un pointeur sur les deux 
triangles adjacents, et les deux sommets qui forment 1' arete, 

- pour chaque triangle, on determine les pointeurs sur les 
trois sommets qui forment le triangle, les trois aretes du 
triangle, les trois triangles adjacents a chaque arete, ainsi 

30 que trois valeurs donnant l'indice d'adjacence du triangle 
courant pour chaque triangle adjacent. 

L' orientation des trois sommets dans le triangle se fait comme 
suit • 

35 

Pour que le maillage global soit correctement oriente, il faut 
que chaque arete soit parcourue dans deux sens differents 
lorsque l'on parcoure chacun de ses deux triangles adjacents. 
En se referent St la figure 4A, 1' arete Ai d'indice i est 
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1' arete opposee au sommet Si d'indice i dans ce triangle. Le 
triangle adjacent d'indice i est le triangle voisin TV, 
adjacent a 1' arete Ai. Enfin, chaque triangle a toujours trois 
triangles adjacents. Pour representer des surfaces a bords (un 
plan a un bord, un cylindre a deux bords) on utilise la notion 
de sommets et de triangles dits « virtuels » (figure 4B) . "Les 
triangles virtuels ne sont pas affiches et la position du 
sommet virtuel n'a pas de consequence sur la geom6trie du 
maillage. 

On donne ci-apres un exemple, en langage C++, d'une attribu- 
tion d'indice. 



Si on appelle neighbor [i], le triangle adjacent numero i, 
15 alors l'indice d'adjacence « whichNeighbor[ i] » est definie 
par une instruction telle que : 

Neighbor[i]->neighbor[whichNeighbor[i] ]=this ; 

Pour des motifs de maillages t^traedriques , la structure de 
20 donn^e d'un volume de t^traedre consiste en une liste de 
t^traedres, une liste de triangles, une liste d'arStes et une 
liste de sommets . 

Pour chaque sommet, on stocke sa position, un tetraedre 
25 adjacent. 

Pour chaque arete, on stocke les pointeurs sur les deux 
sommets adjacents et un pointeur sur un triangle adjacent a 
cette arete. 

30 

Pour chaque triangle, on stocke les pointeurs sur les trois 
sommets qui forment le triangle et les pointeurs sur les deux 
• t^traedres qui sont adjacents au triangle. 

35 Pour chaque t^traedre, on stocke : 

- les pointeurs sur les quatre sommets qui forment le 
t^traedre, 

- les pointeurs sur les six ar€tes du t6trafedre, 

- les pointeurs sur les quatre triangles du t6traedre, et 



20 
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- les pointeurs sur les quatre tetraedres voisins opposes a 
chacun des quatre sommets. 

Chaque tetraedre est oriente. A partir de 1 ' orientation des 
5 quatre sommets ABCD , on cherche 1 ' orientation des quatre 
triangles : BCD, ADC , ABD, ACB. Pour que le maillage soit 
correctement oriente, il faut que chaque triangle soit 
parcouru dans deux sens dif ferents en parcourant chacun de ses 
deux tetraedres adjacents. Enfin chaque tetraedre est tou jours 
10 adjacents a quatre tetraedres. On introduit la notion de 
sommets dits «virtuels» et de « tetraedres virtuels» pour 
prendre en compte les bords du maillage. 

Pour simplifier la structure de l'objet deforme, on procede 
15 a une decoupe de maillage en enlevant de la matiere (plus 
precisement un sommet, en imposant neanmoins que ce sommet ne 
peut pas etre un sommet virtuel). 



Dans le cas d'un maillage triangulaire, si ce sommet n'est pas 
adjacent a un sommet virtuel, alors on considere ce sommet 
comme etant virtuel ainsi que tous ces triangles adjacents. 
Si le sommet est deja adjacent a un sommet virtuel, alors on 
enleve 1' arete qui relie ce point au sommet virtuel, et on 
force tous les triangles adjacents au sommet virtuel a etre 
25 virtuel (figures 5A et 5B) . 

Avec un motif de maillage tetra6drique, on procede en enlevant 
un ensemble connexe de tetraedres. Une fois les tetraedres 
enleves, le maillage doit rester conforme, comme decrit ci- 
apres. Pour cela, il peut etre avantageux d'enlever plus de 
tetraedres que ceux presents dans 1' ensemble de tetraedres 
indique au depart. On peut aussi ajouter ou enlever des 
sommets virtuels lorsque l'on cree des composantes connexes 
ou que l'on supprime des bords. 



30 



35 



Le Demandeur s'est confronte a la difficult^ d' assurer une 
conformite du maillage, dont on donne une definition ci-apres. 
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Pour qu'un maillage triangulaire soit conforme, chaque arete 
doit appartenir au plus a deux triangles et le voisinage 
topologique de chaque sommet doit rester homeomorphe a un 
cercle (ou a un demi-cercle) , c'est-a-dire que, si deux 
triangles partagent un meme sommet, alors il doit exister un 
ensemble de triangles adjacents entre eux qui permettent de 
passer d'un triangle a un autre. En se referent a la figure 
6, le maillage triangulaire represent^ est non conforme car 
les deux triangles ont une intersection suivant un sommet, 
tandis qu'il n'existe pas de triangles adjacents permettant 
de passer de l'un a 1' autre de ces deux triangles. 



Pour qu'un maillage tetraedrique soit conforme, il faut 
que chaque triangle appartienne au plus a deux tetraedres et 
L5 que le voisinage topologique de chaque sommet soit homeomorphe 
a une sphere (ou une demi-sphere s'il est sur le bord).. Il 
faut done que, si deux tetraedres partagent un meme sommet, 
alors il existe un ensemble de tetraedres adjacents entre eux 
qui permettent de passer d'un tetraedre a un autre. 

20 

L* equation qui gouverne le mouvement du modele elastique 
lineaire est liquation dif f erentielle Newtonienne El de 
1' annexe. Cette equation est reliee aux Equations differen- 
tielles rencontrees en mecanique des milieux continus 
25 (equation E2) . 

D'apres la theorie des elements finis, la matrice de masse M 
et la matrice d' amortissement C sont creuses et dependent des 
proprietes de chaque tetraedre. De plus, ces matrices M et C 

30 peuvent etre considerees ici comme diagonales, c'est-a-dire 
que la masse et les effets d* amortissement sont concentres sur 
les sommets du maillage. Cette simplification, appelee masse- 
lumping, decouple les mouvements de chacun des noeuds et 
permet done d'ecrire 1' equation precedente E2 comme un 

35 ensemble d' equations dif f erentielles independantes du type El, 
avec une equation pour chaque sommet. 



Le choix d'un schema d * integration explicite mene a exprimer 
les forces eiastiques a 1* instant t, pour calculer la nouvelle 
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position du sommet a 1' instant t+1, avec liquation E3 de 
1' annexe. Avantageusement , cette equation E3 est valable pour 
tout type de force, lineaire ou non. 

Les schemas explicites sont seulement conditionnellement 
stables, ils convergent done moins vite que les schemas 
implicites . 

Une autre difficulty s'est presentee : pour que le schema soit 
stable, le pas de temps utilise pour 1 • integration numerique 
doit etre suffisamment petit. Le pas de temps critique depend 
de la plus grande valeur propre de la matrice de rigidite et 
des valeurs locales de la masse et de 1 ' amortissement . Pour 
optimiser le pas de temps, la matrice de masse est fixee comme 
15 etant egale a l'identite, et 1' optimisation des valeurs 
locales de 1 ' amortissement , afin de limiter les risques 
d' oscillations, ont pour effet d'ameliorer le conditionnement 
du systeme lineaire, et done d'en stabiliser la resolution. 

20 On se refere maintenant a la figure 7, sur laquelle un 
organigramme represente schematiquement les dif ferentes etapes 
de traitement appliquees a une surface tridimensionnelle 
triangulee, a titre d'exemple. 

25 Le traitement commence par une premiere etape 10, en laquelle 
on applique la structure des donnees decrites ci-avant, avec 
une attribution d* indices adequats. On considere ensuite 
chaque triangle T n a 1' etape 12, pour lequel on cherche les 
coefficients de Lame X n et u„, ainsi que l'aire du triangle A 
30 (etape 14). Pour chaque arete du triangle (j etant un entier 
de 0 a 2), on calcule la difference des carres des longueurs 
d' aretes, entre la longueur courante lj et la longueur d' arete 
precedente L 3 (etape 16). De ces differences de longueurs 
d' aretes au carre, on deduit les forces exercees en chaque 
35 sommet S 3 du triangle T n (force Fj n calculee a 1' etape 18 de 
1' organigramme) . En passant aux triangles suivants, on integre 
les forces sur tous les triangles en bouclant sur l'indice n 
(etape 20) . 
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Aux stapes 30 et 32, on determine les nouvelles positions Q' n 
des sommets S t des triangles deformes, a partir de leur 
position precedente Q n et des forces F n . En bouclant sur 
l'indice n a 1' etape 36, on realise une integration sur tous 
5 les ob jets deformes. L' estimation des nouvelles positions des 
sommets S jr s'effectue preferentiellement en utilisant les 
equations El, E2 et E3 pr#cit6es. 

Dans 1 • organigramme de la figure 7, on a represent^ les deux 
.0 integrations respectives (des forces Ff et des positions des 
sommets S, dans les objets deformes), separement. Bien 
entendu, il peut etre prevu de realiser certaines des etapes 
12 a 20, en meme temps ou apres certaines des etapes 30 a 36. 
Ce choix depend de la structure des deux modules decrits ci- 
15 avant, que comporte le calculateur precite, et, le cas 
echeant, de la structure des donnees. 

II est preferentiellement prevu un test 42 pour savoir si 
1' ob jet a entierement retrouve sa forme de repos, a savoir : 
20 existe-t-il une arete dans l'objet, dont la longueur est 
sensiblement differente de celle de la longueur precedente ? 

En parcourant les triangles, si le calculateur ne trouve pas 
une telle arete, c'est que l'objet est revenu a sa forme de 
25 repos initiale et le traitement s'arrete en 40. 

Sinon, une image de l'objet est affichee par le moniteur VISU 
precite et une boucle de temps (etape 44) est incrementee , 
pour passer a un temps suivant et montrer un ret our progressif 
30 a une forme de repos de l'objet, en utilisant, le cas echeant, 
liquation E3 . Cette boucle de temps peut dtre initialise 
apres la premiere etape 10 de depart. 

Bien entendu, 1 ' invention n'est pas limitee a la forme de 
35 realisation decrite precedemment a titre d'exemple, elle 
s'etend a d'autres variantes d^finies n^anmoins dans le cadre 
des revendications ci-apres. 
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Ainsi, on comprendra que 1' invention ne se limite pas a une 
application a la simulation de la deformation de tissus 
corporels mous. Comme decrit precedemment , le deplacement d'un 
objet rigide (par exemple pour un os solidaire d'un muscle), 
5 solidaire du tissus mous, peut aussi etre simule. II convien- 
dra, le cas echeant, d'utiliser un torseur (d' expression 
connue en soi) , plutot qu'un tenseur de deformation du type 
decrit ci-avant. 



10 La formule de l'energie de deformation W, aussi bien dans une 
maille triangulaire que tetra«§drique, est donnee ci-avant a 
titre d' exemple et depend du materiau dans lequel est prevu 
1' objet dont la deformation est simulee. Neanmoins, elle 
s'exprime de facon generale par une combinaison des traces du 

15 tenseur E et du tenseur E 2 . Dans des cas particuliers, seule 
la trace du tenseur E (ou du tenseur E 2 ) peut intervenir. 
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Re vend ic at ion s 

1. Dispositif de simulation de deformations de materiaux, 
notamment de tissus corporels mous f comprenant : 

- une zone-memoire ( MEM , MT ) apte It stocker, pour un objet au 
moins : 

* des param&tres mecaniques du materiau de 1'objet , 

* des donn^es de position de 1' objet (Lj), prises en des 
sommets (Sj) d'au moins un maillage choisi (Tn), et 

* des donn6es de force (Ff), representant , en intensity et 
position , une sollicitation 2l exercer sur 1'objet, 

- ainsi qu'un calculateur ( \iP , MT ) capable de coop^rer avec la 
zone-m6moire (MEM) pour ^valuer de nouvelles positions des 
sommets (Sj), en fonction de la sollicitation exercee et des 
parametres mecaniques du materiau, 

caracterise en ce que le calculateur comprend: 

- un module ( 12 , 14 , 16 , 18 , 20 ) apte Si determiner iterativement , 
pour chaque maille, l'ecart entre la longueur courante (Ij) 
d'une arete de la maille et sa longueur precedente et/ou au 
repos (Lj), ainsi qu'& stocker des donnees de force respecti- 
ves (Fj n ) pour chaque sommet de la maille, et 

- un module (30,32,36) apte h determiner iterativement , pour 
chaque sommet, de nouvelles donnees de position de ce sommet 
(Sj(Qn)) en fonction de la composition des forces qui s'y 
exerce (Fj n ), de donnees portant sur au moins une position 
anterieure du sommet (Sj(Tn)) et des parametres mecaniques du 
mater iau. 

2. Dispositif selon la revendication 1, caracterise en ce que, 
pour un objet presentant la forme d'une enveloppe tridimen- 
sionnelle creuse, le maillage choisi est triangulaire, de 
sorte que le calculateur est apte ct determiner une composition 
des forces en chaque sommet de triangle, en fonction de 
l'ecart entre la longueur courante de chaque arete du triangle 
et la longueur de cette arete au repos. 

3. Dispositif selon l'une des revendications 1 et 2, caracte- 
rise en ce que, pour un objet presentant une forme tridimen- 
sionnelle pleine, le maillage choisi est tetraedrique, de 
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sorte que le calculateur est apte a estimer une composition 
des forces en chaque sommet de tetraedre, en fonction de 
I'Scart entre la longueur courante de chaque arete du 
tetraedre et la longueur de cette arete au repos. 

4. Dispositif selon l'une des revendications pr6c6dentes, 
caracterise en ce que le calculateur est agence pour determi- 
ner des differences entre les carr<§s des longueurs courante 
(lj) et pr^cedente et/ou au repos (Lj) de chaque arete pour 
determiner ladite composition des forces. 

5. Dispositif selon l'une des revendications precedentes, 
caracterise en ce que la zone-memoire est agencee pour 
stocker r en association avec chaque maille, des paramfetres 
nuScaniques du materiau de la maille (K,V*) , au moins partiel- 
lement definis localement, en particulier au niveau de la 
maille ou d' elements de celle-ci. 

6- Dispositif selon la revendication 5, caracterise en ce que 
le calculateur est agence pour une determiner ledit ecart de 
longueurs courante et au repos, en vue d' estimer une derivee 
de l'energie potent ielle de deformation (IV) de chaque maille , 
exprimee en fonction d'un tenseur de Green-St Venant (E) et 
de coefficients mecaniques propres au materiau dans cette 
maille, 

7. Dispositif selon la revendication 6, caracterise en ce que 
lesdits parametres mecaniques comprennent les coefficients de 
Lame du materiau dans chaque maille consideree (Xn,u n )- 

8* Dispositif selon l'une des revendications precedentes, 
caracterise en ce qu'il comporte un module pour eiaborer une 
structure de donnees et agence pour supprimer des aretes (Ai) 
de maille qui relient deux sommets dits "virtuels" (Si). 

9. Dispositif selon la revendication 8, caracterise en ce que 
le module pour eiaborer la structure de donnees est agence en 
outre pour verifier que le maillage respecte des proprietes 
de conformite predefinies. 
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sorte que le calculateur est apte h estimer une composition 
des forces en chaque sommet de tetraddre, en fonction de 
1'ecart entre la longueur courante de chaque arete du 
tetrafedre et la longueur de cette arete au repos. 

4. Dispositif selon l'une des revendications pr6c6dentes, 
caracterise en ce que le calculateur est agence pour determi- 
ner des differences entre les carr6s des longueurs courante 
(lj) et precedente et/ou au repos (Lj) de chaque arete pour 
determiner ladite composition des forces. 



5. Dispositif selon l'une des revendications pr£c£dentes, 
caracterise en ce que la zone-m§moire est agencee pour 
stocker, en association avec chaque maille / des paramdtres 
15 mecaniques du materiau de la maille (Xn,u n ), au moins partiel- 
lement definis localement, en particulier au niveau de la 
maille ou d ' elements" de" celle^ci . ' ' 



6. Dispositif selon la revendication 5, caracterise en ce que 
le calculateur est agence pour determiner ledit ecart de 
longueurs courante et au repos, en vue d' estimer une derivee 
de l'energie potentielle de deformation (W) de chaque maille r 
exprimee en fonction d'un tenseur de Green-St Venant (E) et 
de coefficients mecaniques propres au materiau dans cette 

25 maille. 

7. Dispositif selon la revendication 6, caracterise en ce que 
lesdits paramfetres mecaniques comprennent les coefficients de 
Lame du materiau dans chaque maille consideree (X^,^). 



8. Dispositif selon l'une des revendications precedentes, 
caracterise en ce qu'il comporte un module pour eiaborer une 
structure de donnees et agence pour supprimer des aretes (Ai) 
de maille qui relient deux sommets dits "virtuels" (Si). 

9. Dispositif selon la revendication 8, caracterise en ce que 
le module pour eiaborer la structure de donnees est agence en 
outre pour verifier que le maillage respecte des proprietes 
de conform! te predefinies. 
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10. Dispositif selon l'une des revendications pr<§c<§dentes , 
caracterise en ce que le module (30,32,36), apte a determiner 
les nouvelles donnees de position des sommets (Q'n) en 
fonction de la composition des forces en chaque sommet, est 
agence pour determiner lesdites nouvelles donnees de position, 
en fonction du temps <E3), ce qui permet de suivre une 
evolution temporelle des positions respectives des sommets. 

11. Dispositif selon la revendication 10, caracterise en ce 
) que la determination des nouvelles donnees de position des 

sommets (Q'n) utilise un modele de resolution d' equation 
differentielle appliquee a chaque sommet. 

12. Dispositif selon l'une des revendications precedentes, 
5 caracterise en ce que le calculateur est capable de determiner 

repetitivement les donnees de position des sommets du 
maillage, en vue de determiner une evolution temporelle (t+At) 
desdites positions. 

0 13. Dispositif selon la revendication 12, caracterise en ce 
qu'il comporte une interface de visualisation (VISU,IG) propre 
a representer l'objet suivant une forme predeterminee et 
maillee selon le motif choisi, et en ce que 1' interface de 
visualisation est apte a cooperer avec le calculateur pour 
25 representer la forme de l'objet en mouvement. 

14. Dispositif selon l'une des revendications precedentes, 
caracterise en ce qu'il comporte une interface utilisateur 
(IU) munie d'un organe de saisie (CLA,MO) pour simuler une ou 
30 plusieurs forces globalement exercees sur l'objet. 
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